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THE UNIVERSAL INITIAL MASS FUNCTION IN THE XUV DISK OF M83 
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ABSTRACT 

We report deep Subaru Ha observations of the extended ultraviolet (XUV) disk of M8J@. These 
new observations enable the first complete census of very young stellar clusters over the entire XUV 
disk. Combining Subaru and GALEX data with a stellar population synthesis model, we find that (1) 
the standard, but stochastically-sampled, initial mass function (IMF) is preferred over the truncated 
IMF, because there are low mass stellar clusters (1O 2_3 M0) that host massive O-type stars; that (2) 
the standard Salpctcr IMF and a simple aging effect explain the counts of FUV-bright and Ha-bright 
clusters with masses > 10 3 M Q ; and that (3) the Ha to FUV flux ratio over the XUV disk supports the 
standard IMF. To reach conclusion (2), we assumed instantaneous cluster formation and a constant 
cluster formation rate over the XUV disk. The Subaru Prime Focus Camera (Suprime-Cam) covers 
a large area even outside the XUV disk - far beyond the detection limit of the HI gas. This enables 
us to statistically separate the stellar clusters in the disk from background contamination. The new 
data, model, and previous spectroscopic studies provide overall consistent results with respect to the 
internal dust extinction (Ay ~ 0.1 mag) and low metallicity (~ 0.2Z Q ) using the dust extinction curve 
of the Small Magellanic Cloud. The minimum cluster mass for avoiding the upper IMF incompleteness 
due to stochastic sampling and the spectral energy distributions of O, B, and A stars are discussed in 
the Appendix. 

Keywords: Galaxies: star clusters: general — Stars: luminosity function, mass function — Galaxies: 
individual:M83 



1. INTRODUCTION 

The initial mass function (IMF) is of central im- 
portance in modeling galaxy formation and evolu- 
tion. Its universality has been supported observatiq nally 
(e.g., lKroupal[2002t iBastian. Covev fcMeverl 120 10ft an d 
adopted in many applications (e.g., iKennicutt] 1 19981 ). 
Observations of resolved stellar populations and inte- 
grated colors of most galaxies are consistent with a uni- 
versal IMF, but a non-universal IMF has been inferred 
recent ly, possibly depending on eny ironmental param- 
eters (|Bastian. Covev fc Meven l20Tot) . The IMF could 
be steeper (or truncated), as compared to the stan- 
dard one, in low density environments, such as dwarfs 
and low surface brightness (LS B) galaxies, thus pro- 
ducing fewer high- mass stars (|Hoversten~fe Glazebrook 
120081; IMeurer et all [20091: iLee et al.ll2009D . On the other 
hand, a top heavy IMF is suggested in dens e starburst 
galaxies (|Rieke et al.lll993t iParra et al.ll2007[ ). 
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Pfla mm-Altenburg fc Krou pa (2008) explained the ap- 
parent non-universality as a purely statistical effect; 
the intrinsic IMF shape is invariant, but the maximum 
mass of the stars in a stellar cluster is limited by the 
mass of the parental cluster, the mass of which is in 
turn controlled by star formation activity. Thereby, the 
IMF integrated over a galaxy - integrated galactic IMF 
(IGIMF ), may apparently dif f er in l ow density environ- 
ments. iKrumholz fc McKeel (|2008l ) suggested another 
model: a low gas density could result in fragmentation 
into gas clumps too small to allow high-mass star forma- 
tion. The non-universality of the IMF potentially alters 
our understanding of star formation and galaxy evolu- 
tion. Its confirmation is urgent. 

The Galaxy Evolution Explorer (GALEX) provides a 
means to investigate the high-mass end of the IMF. The 
combination of GALEX and Ha data can separate the 
populations of O and B stars; the two GALEX bands, 
FUV (1539A) and NUV (2316A), are sensitive to both 
stellar clusters with O stars and with B stars, while Ha 
data are responsive pred ominantly to clusters with O 
stars. IMeurer et aLI (|2009D analyzed the Ha to FUV flux 
ratio over galaxies, found a deficit of ionizing O stars, 
and concluded that the IMF is steeper (or truncated) in 
dwarfs and LSBs. This analysis, however, is sensitive 
to the procedure of dust extinction correction (see alter- 
native analysis; iBoselli et alj|2009l ). Another approach, 
the comparison of Ha flux with clus ter mass, points to a 
universal IMF (|Calzetti et al.ll2010l ). 

GALEX finds UV emission far beyond the optical 
edge of gala xies, i.e., their extended ultraviolet disks 
(XUV disks; IThilker et al1l2005l : iGiTde Paz et al.11200.1 
iThilker et all 120071) . Star formation in such low den- 
sity outskirts is a prime target to study the IMF vari- 
ation. In fact, the radial UV profile extends beyond 
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the Ha truncation radius (Boi ssier et al.l 120071 ) . indicat- 
ing that O stars are extremely rare while B stars are 
abundant. The absence may be an indication of a non- 
universal IM F, or could simply be the effect of mixed 
stellar ages (|Thilker et al.l 120051) - as O stars are the 
primary sources of HII regions and have much shorter 
lives (~ 6 Myr) than B stars (~ 1 00 M yr). Indeed, 
IGoddard, Kennicutt fc Ryan- We ber (2010) averaged Ha 
and FUV fluxes over XUV disks and found the ratio 
consistent with the prediction from the standard IMF, 
although the sensitivity and area coverage of their Ha 
data were limited. 

We present Ha observations of M83 obtained 
with the Subaru Pr ime-focus Camera (Supr ime-Cam; 
Mivazaki et al.ll2002D on the Subaru telescope (|Ive et al.l 
2004) . With its high sensitivity and wide field-of-view 
(34' x 27'), we can cover the entire XUV disk and sur- 
rounding area with two pointings (Figure [T]) and can de- 
tect HII regions around individual massive stars. M83 is 
one of the nearest galaxies with a pro minent XUV disk 
(d — 4.5 Mpc or m — M = 28.27 mag; iThim et alJ l2003) 
and is best suited for the detection of individual HII re- 
gions. This study provides the first complete census of 
the UV-bright cluster population with and without HII 
regions. We count the number of blue clusters in the UV 
and bright clusters in Ha and translate the number ra- 
tio into the ratio of durations when a cluster is UV-blue 
and Ha- bright. We compare the ratio with the predic- 
tion from a stellar population synthesis model using the 
standard IMF. 

2. DATA 
2.1. Subaru 

Observatio ns were made in the Rg and NA659 
(Ha) bands (lOkamura e t al J 12001 F Yoshida et ail l200l 
lHavashino et al.l 120031 ) on September 7, 2010. Two 
Suprime-Cam fields cover the entire XUV disk of 
M83, i.e., the whole area with an HI surface density 
above 1.5 x 10 20 cm~ 2 (roughly the detection limit in 
IMiller. Bregman & Wakkcr 2009) and a large surround- 
ing area. Figure [T] shows the area coverage and the 
HI contour (red) on a GALEX FUV image. The two 
fields are marked with black overlapping solid boxes. 
The exposure time for each field is 5x12 minutes and 
5x4 minutes in the NA659 and Rc bands, respectively. 
The yellow contour indicates the edge of the optical disk 
(roughly, the 25 mag/arcsec 2 isophote). The typical see- 
ing was ~ 1". 

Data were reduced in a standard way using the 
Suprime-Cam reduction software (SDFRED; Yag i et all 
120021 lOuchi et al.l 120041 ). We performed overscan sub- 
traction, flat-fielding, distortion correction, background 
subtraction, and mosaicking. The flux calibration 
was performe d against standard field stars (Okc 1990; 
lLandoltl[l9 92). Rc a nd NA659 bands are p ractically at 
the same wavelength (jHavashino et al.|[2003|) . and we do 
not expect any spatial variation in the Rc /NA659 ratios 
of typical foreground field stars. However, there remain 
residual errors at large spatial scales in the ratios (a few 
perent) which were likely introduced in flat fielding. We 
made a correction by making a NA659/Rc ratio map us- 
ing the average ratio of field stars in each 200" x 200" 
grid and by normalizing the NA659 image. The la pho- 




Figure 1. Subaru coverage and definition of regions in M83. 
The background is a GALEX FUV image downloaded from 
Galexview. Solid overlapping rectangles are two pointings with 
the Subaru Prime Focus Camera (Suprime-Cam; the field-of-view 
is 34' X 27', or 44.5 X 35.3 kpc 2 , in one pointing). The inner (yellow) 
contour is roughly at the traditional edge of the optical disk. The 
outer (red) counter is at an HI surface density of 1.5 X 10 20 cm~ 2 , 
using the HI data from IMiller. B regman fc WakkeJ 1 12009 ). We 
call FUV-bright objects inside the HI contour IN objects (exclud- 
ing those in the inner disk - yellow contour), and those outside 
OUT objects. The Ho-bright clusters (NA659 - R c < -1 mag; 
£|4l are also plotted: green for clusters with > 10 3 Mq and cyan for 
the ones with < 10 3 M Q . 

tometric limits with the ~ 5" aperture (i.e., GALEX 
resolution) are 24.90 and 25.56 AB mag in the NA659 
(Ha) and Rc, respectively. 

2.2. GALEX 

We use archival GALEX far -ultraviolet (FUV ) and 
near- ultraviolet (NUV) images (|Bigiel et all 120101 ) . The 
data are obtained from the GALEX data archive through 
Galexview in the multimission archive at STScI 
(MAST). The G ALEX observations and data reduction 
were discussed in iBigiel et all (|2010t ). The FWHM sizes 
of the point spread func tions (PSF) are measured with 
the IDL software ATV (|Barthl I200l1 ) and are 4.7" for 
FUV and 5.2" for NUV. The la photometric limits with 
the ~ 5" aperture are 26.34 and 26.20 AB mag in the 
FUV and NUV-bands, respectively. 

3. METHODS 

3.1. Stellar Population Synthesis Model 

Our analysis is guided by the stellar populatio n syn- 
thesis model STARBURST99 (|Leitherer et al 1999). The 
Padova stellar evolution tracks with AGB stars are 
adopted. Most UV-bright objects in the XUV disk 
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of M83 are likely ste llar clusters (jThilker et al.1 120051: 
IGU de Paz et all 12007). Therefore, we adopt an instan- 
taneous burst approximation to model their photometric 
evolution. 

We will show that the standard Salpeter IMF (Salpeter 
I1955D explains the observations well. Figure [2] demon- 
strates the differences between the standard and trun- 
cated IMFs. The cluster mass is set to IOOOMq for the 
plot, which is roughly the minimum mass of a cluster to 
fully populate the IMF up to an O star (see Appendix^}. 
Note that a factor of 10 increase in cluster mass shifts all 
the magnitude curves in Figure [2] (top) systematically 
upward by 2.5 magnitudes, but it does not change the 
color curves. For these models, we assume the Salpeter 
IMF with a lower cut-off mass of mi — 0.1M©. The 
upper cut-off mass is set to m u — 100, 20, and 3 M , 
which corresponds to the masses o f the most massive O, 
B, and A-type stars, respectively (jCoxl 12000). In other 
words, the last two models represent the truncated IMF 
(or, almost equivalently, the steeper IMF) without O or 
OB stars. We adopt 20% of the solar m etallicity (0.2Zm) 
from the spectroscopic m easurements (jGil de Paz et al.1 
l2007t iBresolin et al.ll2009f >. For a reference, the spectra 
of individual O, B and A stars are discussed in Appendix 

H 

There are quantitative differences between clusters 
with and without OB stars, i.e., between the standard 
and truncated (steeper) IMFs. Clusters without OB 
stars are fainter by 6 mag in FUV for the first 10 Myr 
and never become as blue as FUV-NUV<0.2— 0.3 mag. 
Therefore, blue clusters (<0.2— 0.3 mag) must have O 
and/or B stars. The FUV-NUV color is insensitive to 
the difference between clusters with only B stars (up to 
m u = 2OM ) and with both O and B stars (100M o ). 
This is primarily because the GALEX bands trace only 
the Rayleigh- Jeans side of stellar thermal radiation for 
both O and early- type B stars (Appendix [B]) . The FUV- 
NUV color is therefore a useful indicator of the presence 
of O and/or B stars in a cluster. 

The NA659 — Rc color differentiates clusters with O 
stars from those without them, since O stars emit a 
significa ntly greater number of ionizing photon s than 
B stars ([Sternberg. Hoffmann fc Pauldrachl 120031) . This 
is demonstrated in Figure [5] (bottom). Clusters with 
NA659 — Rc < -1 mag must have O stars, and those 
with NA659 — Rc < -0.3 mag should have O and/or 
B stars. For this plot, we calculated the Ha line lu- 
minosity Lhq using the Lyman photon flux Qi, yc from 
STARBURST99 and the relation L H a[ergs _1 ] = 1.37 x 
10~ 12 <3Lyc[s -1 ] from case B recombination with an elec- 
tron temperature of 10 4 K (|Gavazzi eTaD[2Q03). Case 
B recombination is perhaps a reasonable assumption 
even in the low-density environment, since the optical 
depths in t he Lyman resonance l i nes a re very likely 
quite large (jOsterbrock fc Ferlandl 120061 ) in any star- 
forming conditions. The NA659 magnitude is based 
on the stellar continuum luminosity plus Lh q - We 
ignore the effect of [Nil] emission since the effect is 
very small dGil de Paz et al.ll2007t iKennicutt et alj|2008t 
iGoddard. Kennicutt fc Rvan-Weberll2010L see qCT ]) . 



3.2. Cluster Counts 



t 

I s 



lOOMsun 



20Msun 



I.OZsun 
0.4Zsun 
(UZsun 



mu=3Msun 



Tlrfie [^4af j 



3Msuri 



20Msun 




"IGOMsun 



Hag Timfl \f«Ql) 





JMiun 






20Msun 












lOOMsun 


J 

Jf 











103 Urn. |j»Kirj 

Figure 2. Photometric evolution of stellar clusters from Star- 
burst99. A single-burst model with the Salpeter IMF is adopted. 
A low metallicity (0.2 Zq) is adopted (red) for the analysis, but 
models of higher metalicities (O.4Z0 and I.OZq) are also plotted 
(green and cyan, respectively). The lower mass boundary for the 
IMF is O.IMq, and the upper mass boundary is IOOM0 (including 
all types of stars), 20Mq (no O stars), and 3Mq (no OB stars). 
The adopted distance of M83 is d = 4.5 Mpc. Top: FUV mag- 
nitude for a cluster mass of 1000M Q . Middle: FUV-NUV color. 
Bottom: NA659-i? c color. 

Counting blue clusters in FUV-NUV and NA659 - R c 
constrains the IMF at its high mass end. Clusters be- 
come bright in Ha ("blue" in the color NA659 — Rc) 
predominantly due to O stars, while FUV-NUV is due 
to both O and B stars. Since O stars have much shorter 
lives than B stars, NA659 — Rc changes faster than FUV- 
NUV (Figure [2]). Under the assumptions of instanta- 
neous cluster formation and a constant cluster forma- 
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tion rate, the numbers of blue clusters in FUV-NUV and 
NA659 — Rc (N\rv an d Nn a ) should be proportional to 
the durations that clusters are blue (tuv and ^Ha, re- 
spectively). Thus, we expect Nn a /N\jy — tn a /t\jy if the 
difference in the populations is due to a simple aging ef- 
fect starting from the initial populations of the standard 
IMF. The assumption of a constant cluster formation 
rate should be approximately accurate when integrated 
over the large XUV disk and over a short timescale, i.e., 
the young ages of clusters which can be detected with 
our color selection criteria (below). In fact, the clusters 
that we count in §5.31 are distributed over the UV disk 
(Figure P. 

For the standard IMF (m u = 1OOM and 0.2Z o ), a 
cluster is blue (e.g., FUV-NUV< 0.0 mag) for t uv = 70.8 
Myr and Ha-bright (e.g., NA659 - R c < -1.0 mag) 
for 5.8 Myr. Therefore, the number ratio should be 
Nua/Nxjy = 0.08 if each population is the consequence 
of the cluster aging effect. We will compare the observa- 
tions with this theoretical prediction. 

We identify and count stellar clusters at the GALEX 
resolution (~ 5"; ~ 109 pc at the distance of M83). It is 
possible that clusters may be blended at this resolution; 
however, even in that case, our analysis is valid if clusters 
formed in such proximity are typically coeval. 

4. SAMPLE SELECTION 
4.1. FUV-Selected Objects 

UV-bright objects are identified bas ed on the FUV im- 
age u sing the S Extractor package (|Bertin fc Arnoutsl 
1996). Saturated stars and image edges are masked be- 
fore the identification. Objects within the traditional op- 
tical disk (yellow contour in Figure Q} are excluded. We 
confirmed the very low probability of false detection us- 
ing the negative image technique. SExtractor is run 
against the sign-flipped FUV image and detected only 
one false target at 25.41 mag in FUV. Thus, false detec- 
tions have a negligible impact on our number counts. 

The aperture photometry is carried out in the four 
bands with the dual image mode of SExtractor us- 
ing the FUV image as the reference. The la photomet- 
ric limits with the aperture of the GALEX resolution 
(~ 5") are 24.90, 25.56, 26.34, and 26.20 mag in the 
NA659 (Ha), R c , FUV, and NUV-bands, respectively. 
The sensitivity in NA659 (FWHM=120A) corresponds 
to a flux of 7.98 x 10 34 ergs _1 at a distance d — 4.5 Mpc, 
which is sufficient to detect an HII region around a single 
B star (e.g., 3.44 x 10 35 ergs _1 in case of B0; assuming 
case B recombination; see Appendix [Bj . We note that 
the NA659 band includes stellar continuum emission as 
well as Ha emission. 

The contamination from [NII]A6548, 6583 in the 
NA659 photometry is neglected. The line ra- 
tio [Nil] /Ha is small (< 0.1) in the outskirts of 
galaxies dGil de Paz et alTl2007t iKennicutt et all 120081: 
iGoddard. Kennicutt fc Ryan- Weberl 120101 ) . Even in the 
extreme case that the stellar continuum has zero contri- 
bution to the NA659 flux, the error is only 0.1 mag (and 
likely less in more realistic conditions). 

Figure [1] shows that the Subaru data cover both the 
inside and outside of the HI contour (roughly the HI de- 
tection limit). Recently-formed FUV-bright objects are 
more likely to reside inside the HI contour as star for- 



mation occurs in gas. For the background subtraction 
in §4.41 the objects inside the HI contour (hereafter, IN 
objects) and outside (OUT objects) are separately cata- 
loged. Note that we intensively used the IMCAT software 
package (Nick Kaiser, private communication) for the 
following data analysis. 

In M5 . 31 we will discuss a selection of clusters based on 
their NA659 — Rc color; but in advance, Figure [T] shows 
the locations of clusters with NA659 — Rc < — 1 mag 
(i.e., bright in Ha). They are distributed over the entire 
XUV disk, and all are within the HI contour (see §4.4j) . 

4.2. Extinction Correction 

The extinction curves of [Pel (|1992ft are adopted to 
correct for Galactic and internal extinctions. The ex- 
tinctions in FUV, NUV, R c , and NA659 are (A FUV , 
^nuv, A Rc , A NA659 ) = (2.58 , 2.83 , 0.81, 0.81)A V for the 
Galaxy, (3.37, 2.61, 0.80, 0.80)A V for the Large Magel- 
lanic Cloud (LMC), and (4.18, 2.59, 0.79, 0.79)A V for 
the Small Magellanic Cloud (SMC), where Ay is the ex- 
tinction in V-magnitude. We note that the FUV-NUV 
color becomes bluer under the Galactic extinction, but 
redder under the LMC and SMC extinction. 

The Galactic extinction cor rection is made for a ll ob- 
jects using Ay = 0.218 mag (jSchlegel et al.lll998| ). We 
do not apply an internal extinction correction until |J51 
but it will be corrected using Ay = 0.1 mag and the 
SMC extinction curve for the following reason. The mea- 
sured metallicity of mo st outer objects is low (~ 0.2Z Q ; 
IGil de Paz et all 120071 ). and therefore, the extinction 
curve of the SMC is perha ps the most appropriate . Us- 
ing E(B - V) values from IGil de Paz et all (|2007l) . we 
estimate the amount of internal extinction to be around 
Ay ~ 0.1 mag with outliers (e.g., 0.0 - 1.0 mag). We 
therefore adopt Ay ~ 0.1 mag for the statistical correc- 
tion, and this degree of extinction w ill b e supported by 
a comparison with a model color in ^ 15. 21 

4.3. Completeness 

The mass detection limit is determined in color mag- 
nitude space with the help of models. Figure [3] shows 
color-magnitude diagrams for IN objects. The overlaid 
solid lines are the 3a and 5a magnitude detection lim- 
its and model evolutionary sequences for clusters with 
masses of 10 2 , 10 3 , and 10 4 M Q . The models are for clus- 
ters with the standard Salpeter IMF with (mi, m u ) — 
(0.1, 100M Q ) and metallicity O.2Z . Models with lower 
m u do not reproduce clusters as blue as observed (Figure 
[2 and see discussion below) . Clearly, the mass detection 
limit changes as a function of color in Figure [3l it de- 
creases as clusters become redder. We note that the Rc 
and NA659 bands are practically at the same wavelength, 
and therefore, NA659 — Rc is zero without Ha emission. 
A negative NA659 — Rc color indicates the presence of 
associated Ha emission. 

The top panel shows that the Subaru observations de- 
tect clusters down to ~ 10 2 M Q . Clusters become red- 
der and fainter as they age, and those around the de- 
tection boundary (10 2 Mq) become too faint when they 
become redder than NA659 — Rc ~ —0.6 mag. In other 
words, the sample is complete down to 1O 2 M in the 



See http: / /www. ifa.hawaii.edu/~kaiser/imcat_ 
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Figure 3. Color-magnitude dia gram s of all the objects inside the 
HI contour (IN objects defined in SgjJ. Top: NA659-R C vs NA659. 
5<r and 3cr detection limits are illustrated with solid straight lines. 
The model evolutionary sequences with the standard Salpeter IMF 
are displayed with solid curves for clusters with 10 2 , 10 , 1O 4 M0. 
Bottom: The same as the top panel, but for FUV-NUV vs FUV. 
Single 03 and BO star models are also plotted in both panels. 

range NA659 — Rc < —0.6 mag. It is complete down to 
10 3 M Q in NA659-R C < 0.0 mag. The extinction is neg- 
ligible, as shown by the extinction vector (for Ay = 0.1 
mag) in the plot. 

The bottom panel shows the corresponding color- 
magnitude diagram in the UV. Again, the mass detection 
limit is a function of the color range (and thus, the clus- 
ter age). The detection limit is ~ 10 3 M Q for the color 
FUV-NUV< 0.0 mag and is even deeper (~ 1O 2 M ) for 
FUV-NUV< —0.2 mag. The extinction is not negligible 
in FUV and NUV and will be discussed later. 

Note that the very blue objects in the bottom panel 
(FUV-NUV< -0.3 mag and FUV - 22 - 24 mag) are 
bluer than any model cluster with the standard IMF can 
reproduce, and are located approximately in the region 
where single O stars appear (the models of 03 and B0 
stars are also plotted). These objects are not included 
in our number count, since they appear below the mass 
selection criterion in §5.31 We will discuss these objects 
in g5U 

4.4. Background Subtraction 

The contamination of background (and foreground) 
objects has been an obstacle i n identifying the 
objects within nearby galaxies (iWerk et al.l 120101 : 
iGoddard. Kennicutt fc Ryan- Weber] 12010ft . The large 
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Figure 4. NA659-Rc vs FUV-NUV color-color plot of the objects 
above 5cr-detection and above 10 3 Mq. Lines are model predictions 
with the Salpeter IMF with m u = 100M Q for 0.2 Z Q . A reddening 
vector of Ay = 0.1 mag from the SMC extinction curve is also 
plotted. We assume that stellar clusters (traced in Rc) and the 
surrounding ionized gas (in NA659) suffer from extinction by the 
same amount. Top: FUV-bright objects within the HI contour (IN 
objects). Bottom: the FUV-bright objects outside (OUT objects). 

field-of-view of Suprime-Cam is a clear advantage as it 
enables us to remove such contamination statistically. 
The objects outside the HI gas contour (OUT objects; see 
H4.1[) are most likely in the background, whereas IN ob- 
jects include stellar clusters within M83 as well as back- 
ground objects. 

The IN and OUT objects show different colors statis- 
tically. Figure H shows the plots of NA659-R C vs FUV- 
NUV for the IN and OUT objects above the 5a mag- 
nitude detection limit and above the 1O 3 M0 threshold 
line in Figure [3] (bottom) . Figure [5] is a histogram of 
FUV-NUV color for the IN (green) and OUT objects 
(blue). Clearly, the IN objects are bluer in FUV-NUV 
than the OUT objects and preferentially have FUV- 
NUV<0.2-0.3 mag, while the OUT objects are more 
likely to have FUV-NUV>0.2-0.3 mag. Figure H shows 
that virtually all objects with NA659 — Rc < -0.3 mag 
(OB stars) are IN objects; almost no young stellar clus- 
ters exist outside of the HI gas. This clear distinction 
ensures that the objects outside of the HI gas are not 
the young clusters in M83 and represent the background 
population. Note that IN objects are young stellar clus- 
ters, not very old clusters with many evolved FUV- 
bright stars, such as extreme horizontal branch stars, 
blue stragglers or planetary nebulae. The old clusters 
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Figure 5. Histograms of FUV-NUV color of the objects above 5a- 
detcction and above 10 3 Mq. Colors indicate the objects within the 
HI gas (green), the ones outside of the HI gas (background objects; 
blue), and the ones associated with M83 (red). The associated ones 
are calculated by subtracting OFF from ON with scaling for the 
area coverage difference. 

do not become as blu e in FUV-NUV as our IN objects 
(jDonovan et al.ll2009D . 

The background contamination in the IN objects can 
be removed statistically; we subtract the OUT objects 
from the IN. In Figure [5] the number of the OUT ob- 
jects are scaled by the IN/OUT area ratio (55%) and 
subtracted from that of the IN objects. The derived his- 
togram (red) represents the population of young clusters 
with mass > 10 3 M Q in the XUV disk. 

Virtually all clusters after the subtraction show blue 
color in FUV-NUV (< 0.2-0.3 mag; Figure [5]), suggest- 
ing that all of these young clusters have O or, at least, B 
stars (see Figure [2]). Thus, massive stars are populated 
even in the diffuse gas environment. It is also remark- 
able that almost no red objects (FUV-NUV >0. 2-0.3 
mag) exist. This confirms that OUT objects are dis- 
tributed everywhere across our ficld-of-view and repre- 
sent the background objects; thus, our background sub- 
traction is fairly successful. 

There is a possibility that old stellar clusters (i.e., 
FUV-NUV >0.2— 0.3 mag before extinction correction) 
exist both in the IN and OUT regions and are removed 
from the IN population. However, it should not affect 
our cluster count analysis as we discuss only the young 
clusters. Such a population, if it exists, must be fairly 
uniformly distributed both in the IN and OUT regions, 
since the count after the subtraction is nearly zero in 
FUV-NUV >0.2-0.3 mag. We also note that the ap- 
parent magnitudes of low mass clusters (< 1O 3 M ) fall 
below our detection limit after ~ 100 Myr (Figure [5] top). 
Again, this is not a problem with the selection criteria 
that are adopted in our analysis ( fc|5.3[) . 

5. DISCUSSION 

5.1. Truncated vs Stochastic IMFs 

The presence of ionizing stars (O stars) is evident even 
in small clusters (< 10 3 M Q ) in Figure[31 suggesting that 
the IMF is not truncated. Massive stars do not form un- 
der the truncated IMF at all. On the other hand, they 
should appear occasionally even in small clusters under 
the stochastic IMF, though they tend not to appear in 
a small sample due to their low probability of forma- 



tion in small clusters. Such a stochastic effect becomes 
particularly apparent when the cluster mass is small 
(< 1000M©). Our Subaru observations detect small clus- 
ters (down to ~ 10 2 M Q ) and confirm that even those 
small clusters show Ha excess (i.e., NA659 — Rc < 0.0 
mag); and some even have NA659 — Rc < — 1.0 mag, a 
color which only O stars can produce. Thus, O stars are 
present in the 10 2_3 M Q clusters. The IMF is not trun- 
cated even in the low-density and low-metallicity envi- 
ronment. 

More precisely, it is necessary to account for the 
stochastic effects in discussions of very low mass clusters 
(< 10 3 M Q ) - and even then, our result still supports the 
stochastic IMF over the truncated IMF, as we discuss 
below. We first need to consider the likelihood that a 
cluster with a given small mass has a massive star. In 
Appendix |A1 we estimate that a 10 3 M Q cluster is likely 
to have as massive star as 40M Q ; the cluster mass must 
be greater than ~ 400M© to have a high probability of 
hosting at least one O star (> 20M Q ); a 10 2 M Q cluster 
likely hosts up to a ~ 1OM star. 

The NA659 — Rc color is sensitive to the mass of the 
most massive star in a cluster when the cluster is young, 
whereas FUV-NUV does not distinguish O stars from 
B stars well ( H'S.ll Figure O see also Appendix IB1 . If 
a cluster has no massive star, it should appear redder 
(NA659 — Rc > -1 mag; Figure^]) than the observed clus- 
ters in the range 10 2_3 M Q ; yet some clusters appear in 
the range NA659 — Rc = -l~ -2 mag. Massive O stars in 
small clusters should be populated stochastically. In ad- 
dition, a stochastically-populated single O star can out- 
shine a cluster with a small mass (1O 2_3 M0; Figure [2]). 
This may result in an overestimation of the cluster mass 
when we apply the standard IMF with m u = 100M e to 
the cluster (as done in Figure [3]). Therefore, the masses 
of small clusters that appear in the range 10 2_3 Mq in 
Figure [3] could be actually even less. This prefers the 
stochastic IMF, as even the less massive clusters have O 
stars. 

The discussion so far is about clusters of stars. An- 
other possible evidence for the presence of O stars is the 
population of very blue objects (FUV-NUV< —0.3 mag) 
in Figure [3] (bottom). They have the apparent magnitude 
and color similar to those of a single O star. They are 
possibly isolated, single O stars or small clusters with 
an O star(s) but without many low-mass stars. Higher 
resolution images are necessary to confirm them; their 
presence, if confirmed, may be stronger evidence for O 
stars that are formed in the low-density environments. 

5.2. Internal Extinction Correction 

The internal extinction is assumed to be Ay = 
0.1 mag from spe ctroscopic measurements ( M.'2l 
IGil de Paz etalll200l . This de gree of extinction on av- 
erage is also supported by a comparison with the model 
curve in Figure [4] (top). The extinction vector, corre- 
sponding to Ay = 0.1, is also plotted. An extinction 
correction by this amount creates agreement between the 
UV data and the model curve, providing additional sup- 
port for the adopted Ay. 

We assume that a stellar cluster (traced in Rc) and the 
surrounding ionized gas (in NA659) suffer from extinc- 
tion by the same amount. The two components may have 
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different spatial distributions and a slightly different at- 
tenuation (e.g., E(B - V) con tiraium ~ 0A4E(B- V) eas in 
case of more active star-forming environments; iCalzettil 
120011 ). However, the difference is negligible for our anal- 
ysis since the extinction is small (Ay = 0.1). 

5.3. Universal IMF 

Detected clusters with a color of NA659 — Rc < — 1 
mag are evidence that some clusters host O stars. Here 
we test whether their abundance, relative to the ones 
without O stars, can be explained by a simple aging ef- 
fect with the standard IMF. Under the assumptions of 
instantaneous formation of star clusters and a constant 
cluster formation rate over the XUV disk, the numbers 
and durations of blue clusters should show the relation 
Nn a /N\jv = tua/t\jv, if the relative cluster populations 
are due to the aging effect. 

Consistent sample selection between GALEX and Sub- 
aru data is crucial in this analysis. Each of the four 
photometric bands has its own detection limit (Figure 
[3]), and a bias could be introduced if only magnitude 
detection limits are taken into account. We therefore 
use the cluster mass and color as thresholds for cluster 
counts. As discussed in §4.31 both GALEX and Subaru 
data are complete at > 5a, in colors NA659 — Rc < —0.3 
mag and FUV-NUV< 0.0 mag down to a cluster mass 
10 3 M©. This is true even after the internal extinction 
correction (£ !5.2j) if a slightly lower detection limit (3.4a) 
is adopted. The extinction correction is made for all four 
bands, which is equivalent to shifting the data points 
along the extinction vectors in Figure [3] We select the 
mass limited sample after the extinction correction using 
the technique discussed in Sj4j 

Almost all of the clusters are in the range FUV-NUV < 
mag when the extinction is taken into account (Figure 
El; thus, we adopt this color range for our analysis. This 
is the range where B stars definitely, and O stars possibly, 
exist in the clusters; we count how many of them have 
O stars below. Clusters are blue, FUV-NUV< mag, 
for only a short duration. With the standard IMF of 
(mt,m u ) = (0.1,100M Q ) and 0.2Z Q , this is i uv = 71 
Myr (Figure [2]). The short t\jy is advantageous, since 
the assumption of a constant cluster formation rate is 
likely valid for a short period over the large area (i.e., 
the entire XUV disk) - no local events likely bias the 
statistics. On the other hand, clusters with O stars show 
NA659 - R c < -1 mag for only 5.8 Myr. The model 
therefore predicts ^hoAuv ~ 0.08. 

The number of blue clusters with FUV-NUV < mag 
is counted using Figure El (red) - after the internal ex- 
tinction correction, it is Atjv = 88 ± 9. From Figure El 
N Ha = 9 ± 3 for NA659 - R c < -1 mag (see Figure EJ). 
The errors are based on the Poisson distribution, and the 
counts are made after the internal extinction correction. 
Therefore, the observed ratio is -/Vh^/AW ~ 0.10 ±0.03. 
This is consistent with the prediction of the standard 
IMF (within a la deviation). The counts of stellar clus- 
ters support the standard IMF with a simple aging effect. 

It is worth noting that the slight deviation does not 
likely have a statistical significance. However, if any- 
thing, it is inclined toward the top-heavy IMF side, as 
opposed to the previous suggestion of a truncated IMF. 

The combination of the parameters we adopt [i.e., 
low metallicity (0.2Z Q ), SMC extinction, and Ay = 0.1 




Figure 6. NA659 images of the Ha-bright clusters with NA659 — 
Rc < — 1 mag and > 10 3 Mq. O-type stars must be present to 
produce this cluster color. The size of each image is 24". 2 X 24". 2. 

mag] is the best to consistently explain the new data, 
model, and previous spectroscopic studies. The result 
for the standard IMF does not change even if a higher 
metallicity (~ O.4Z ), and hence the LMC extinction 
curve, is adopted. If we ignore the spectroscopic stud- 
ies, the comparison of the data and model (e.g., Fig- 
ure 2|) provides Ay ~ 0.2 mag. With this extinction, 
we count iVuv = 86 ± 9 and Nn a = 11 ± 3 and obtain 
A Hq /A uv ~ 0.13±0.04. The model ratio is 0.10 with the 
metallicity 0.4Z Q , where t^ a = 5.2 Myr and t\jy = 51.3 
Myr. Therefore, the observations and model are con- 
sistent even if a different metallicity and extinction are 
assumed. 

5.4. Ua-to-FUV Flux Ratio 

The Ha-to-FUV flux ratio (^Ha//iyFUv) over the en- 
tire XUV disk is another measure of the high-mass end 
of the IMF. We find that the observed ratio is consis- 
tent with the one predicted from the sta ndard IMF. A 
lower ratio may sugge st a truncated IMF (|Meurer et alJ 
l200l iLee etafl [20091 ). th ough an extinction correction 
introduces an ambiguity ([Boselli et al.l [2009) . We can 
estimate the ratio more accurately since we derived the 
extinction of Ay — 0.1 mag statistically ( §5.21) . 

In a Starburst99 simulation with a constant star for- 
mation rate, a metallicity of 0.2Z Q , and the Salpeter IMF 
with (mi, m u ) = (0.1, 100)M©, the ratio approaches an 
approximate constant and is, e.g., log(-FH a //„FUv) ~ 
13.10 in unit of log(Hz) after 1 Gyr assuming a star for- 
mation rate of 1M Q /yr (though the constant does not de- 
pend on these values). IFumagal li. da Silva fc Krumhold 
(|2011h calculate the ratio with their own model using the 
Kroupa IMF and derive ~ 13.22 (from their Figure 1). 
The scatter between the models is about 30%. Compared 
to these values, th e observed rat i os pre viously reported 
in LSB galaxies by lMeurer et al.l ()2009f ) - as low as one 
third of the values calculated in these models - deviate 
at a significant level with respect to the model uncertain- 
ties. 

We calculate the Ha to FUV flux ratio over the XUV 
disk of M83 by summing up the fluxes of all > 3a 
detections in the FUV image f £)4.ip . Assuming that 
OUT objects represent the background population ( H4.4p 
and subtracting their total fluxes form those of IN ob- 
jects (accounting for the IN/OUT area ratio), we ob- 
tain F Ha = 8.36 x lO^ergs" 1 and /„fuv = 7.23 x 
10 25 ergs -1 Hz -1 . Adopting Ay = 0.1 mag and the 
SMC extinction curve ( ij4.2p . the intrinsic ratio after 
the extinction correction is log(FH a //j/FUv)int = 12.93. 
This value is only ~ 30% lower than our model pre- 
diction, and not as much as what the previous study 
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found (|Meurer et al.l 120091) . As the difference between 
this value and our model calculation is similar to the 
scatter between the models, we consider a truncated IMF 
to be unwarranted in the case of M83. 

6. SUMMARY 

Using new deep Subaru Ha data and archival GALEX 
data of the extended ultraviolet (XUV) disk of M83, we 
show that the standard Salpeter IMF explains the num- 
ber counts of stellar clusters in the low density, low metal- 
licity environment. Even in clusters with small mass 
(10 2 ~ 3 M Q ) O stars are populated, at least occasionally, 
suggesting that the stochastic IMF is preferable to the 
truncated IMF. The Ha to FUV flux ratio over the XUV 
disk also supports the standard IMF. 

The analyses were guided by the stellar population 
synthesis model Starburst99. The new data, together 
with the population synthesis model and previous spec- 
troscopic studies, provide overall very consistent results 
at the metallicity (~ 0.2Zq), internal extinction (Ay ~ 
0.1 mag), and colors of the stellar clusters in the XUV 
disk. We counted the numbers of FUV-bright clusters 
(which should have O and/or B stars) and Ha-bright 
clusters (having O stars), whose ratio constrains the 
high-mass end of the IMF. The ratio of number counts 
can be translated to the ratio of their lifetimes with 
Starburst99, under the two assumptions of instanta- 
neous cluster formation and a constant cluster formation 
rate in the XUV disk. As a result, the number counts 
are consistent with the Standard IMF and a simple aging 
effect. 
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XUV disk), enabling us to remove the background con- 
taminations statistically, which had been a difficulty in 
previous studies. We found that virtually all blue clusters 
(FUV-NUV< 0.0 mag after the extinction correction; or 
< 0.16 mag before the correction) are within the HI gas. 
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APPENDIX 

MINIMUM STELLAR CLUSTER MASS TO AVOID STOCHASTICITY 

We consider the threshold mass of a stellar cluster where the incomplete sampling (filling) of the IMF (i.e., stochastic 
effect) becomes important. It depends on the definition of "stochastic sampling" , because in reality an O star with 
the highest mass ~ 100M Q is rarely found (at least in the Milky Way), and even a massive cluster (e.g., 10 6 M Q ) may 
not host any star with 100M Q in a statistical sense. On the contrary, we could define the fully-sampled IMF as the 
one in which at least one O star of any subclass (> 2OM ) is populated. Of course, a cluster is more likely to host 
a 20M© star than a lOOM© sta r, and the threshold mass should be lower in the former case. One could also take a 
different approach; for example, Calzctti et al.l (|2010f ) found ~ 10 3 M Q as the threshold mass where the Ha luminosity 
over cluster mass starts to deviate from a constant. 

Here, we take a simpler approach to illustrate a range of threshold cluster mass that may provide an important 
guideline. The threshold mass necessary to populate the IMF to some upper mass cutoff depends on that maximum 
stellar mass 

The stellar initial mass function (\P(to)) with the upper and lower stellar mass cuts (m u and mi, respectively) is 
normalized as 

Mci = / m^(m)dm, (Al) 

J mi 

where M c \ is the total mass of a stellar cluster 
coefficient (A) is 

A 

If we sum up the masses of all stars with masses greater than a threshold stellar mass m t , the sum is 

pm u 

M(> m t ) = / m^>(m)dm. (A3) 

J rn t 

In order to populate at least one star as massive as mt, this sum should be greater than mt [i.e., M(> mt) > m t ]; and 
therefore, 



Assuming a power-law IMF (>F(m) = Am a ) the normalization 

-M d 2 J~° 2 _a - (A2) 
m u — m, 



Subaru H q Observations of M83 XUV Disk 



9 



io 5 




20 40 60 80 100 

mt [Msun] 

Figure 7. Minimum cluster mass to avoid the stochastic sampling effect of the IMF, as a function of the maximum stellar mass that one 
wants to populate unstochastically. The cluster mass needs to be 10 3 , 10 4 , and a few 1O 5 M0 to populate an O star up to 40, 80, IOOMq, 
respectively. 

in the case of a power-law IMF. This condition is necessary to unstochastically populate the IMF up to m t in the 
cluster, providing a threshold cluster mass to have a star up to mass mt- 

A similar calculation also provides a guideline of the maximum stellar mass that a cluster with M c \ is likely to host 
- the IMF would be filled up to m t if it is populated smoothly from the lowest stellar mass mi toward the highest 
until the total mass reaches the cluster mass M c \. Thus, the maximum stellar mass in the cluster is around m t . 

Figure [7] shows the threshold cluster mass M c \ as a function of the maximum stellar mass mt that should be 
populated. We assumed the Salpeter IMF (a = 2.35), mi = 0.1M Q , and m u = 100M Q . In order to have a 40M© star 
the cluster mass should be at minimum 10 3 M Q , while to have 80M© it needs to be 10 4 M Q . All stars with > 20M© are 
classified as O-type stars; and to have at least one O star the cluster mass should be greater than ~ 400M©. To fill 
the IMF literally up to ~ 100M©, the cluster mass should be on the order of 10 5 Mq. If we adopt 50M© (arbitrarily) 
as a typical mass of O stars, the cluster mass should be > 2OOOM0 to fully populate the IMF up to this typical O star 
mass and avoid the stochastic effect. 

These threshold cluster masses are meaningful only in a statistical sense. A stochastic effect could appear even above 
this threshold mass if a Monte-Carlo simulation is performed. Even if a cluster mass is above the threshold and the 
cluster can host a star with m t , this mass could be split and allocated to multiple stars with smaller masses (whose 
total is mt). On the other hand, a cluster with mass below the threshold could occasionally have a star with mass 
> mt at the expense of many less massive stars. The threshold cluster mass should be taken as a minimum mass to 
avoid the stochastic effect, only in a statistical sense. Despite this limitation, our estimate provides a crude guideline 
for discussion of the stochastic sampling of the IMF. 

UV SPECTRA OF O, B AND A STARS 



The spectra of 03, BO, and AO-type stars from the Kurucz's Atlas9 grid (!Kuruczlll992f) are plotted in Figure El 
The masses of t hese stars are approximately 87.6, 17.5, and 2.9M©, respectively ([Sternberg. Hoffmann fc Pauldrachl 
120031: lCoxll2000D . We adopt a metallicity of [Fe/H] = -0.5 (or - 0.3Z©) which is available in the grid. The effective 
temperature (T e ff ) and gravity (log g) are required to find a model from the grid, and the stellar radius R is necessary to 
calculate the total flux. From the Atlas9 grid, we choose the mode ls whose para meter sets are the closest to the stellar 
parameters presented in lSternberg. Hoffmann fc Pauldrachl (pOOl or lCoxl (f2000l ): (T off , \ogg, R) = (50000 K, 5.0, 13.2 
R Q ), (30000, 4.0, 7.4), and (10000, 4.0, 2.4) for 03, B0, and AO, respectively. We note that the spectra of high-mass 
stars are not calibrated very well as only few of those exist in the solar neighborh ood for calibration. In addition, an 
adjustment of the paramete rs is necessary to fit observed spectra of individual stars ([Sternberg. Hoffmann fe Pauldrachl 
2003; iLeitherer et al.l [2010). The spectra that we present here are only a reference to understand our observations. 

The continuum slopes are more or less the same between the O and B stars in the wavelength range of GALEX 
bands. On the other hand, the A star spectrum is rising from FUV to NUV. The FUV absolute magnitude (AB) 
and FUV-NUV color of the individual stars for three metallicity [Fe/H]=+0.0 (1.0Z Q ), -0.5 (0.3Z Q ), -1.0 (O.1Z ) are 
in Table [T] We also list the Lyman continuum photon emission rate (QLyc), which can be used to calculate an Ha 
luminosity (Ljjalergs -1 ] = 1.37 x 10 -12 Qlvc[s -1 ] in case B recombination). The Ql yc is calculated by integrating the 
model spectra over the frequen cy above the ionizing frequency (i.e . , the wavelength of 912A). These Qhyc values are 
4-5 times smaller than those in [Sternberg, Hoffm ann fc Pauldrachl ([2003[ ) , reflecting the uncertainty in model spectra 
of high-mass stars. 
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Figure 8. Top: UV spectra of 03 (87.6Mq), BO (17.5M0) and AO (2.9M0) stars. The y-axis is an apparent magnitude at the distance 
of M83 (d = 4.5 Mpc; m - M=28.27 mag). Bottom: transmission curves of GALEX FUV and NUV bands. 



Table 1 

UV Color of High-mass Star Models 



[Fe/H] = +0.0 -05 -1.0 



Type logQ Lyc FUV FUV-NUV FUV FUV-NUV FUV FUV-NUV 
(s _1 ) (mag) (mag) (mag) (mag) (mag) (mag) 
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Note. — FUV magnitudes are the absolute magnitude. 
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